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ABSTRACT: 

A sum of normal acceleration and pitch rate appears to be the best vari- 
able to use to control aircraft in the longitudinal axis. The C*-Criterion 
specifies that the time response of this quantity called C* must fall in a 
prescribed envelope for all speeds and altitudes. It is equivalent to re- 
quiring the control system to hold the coefficients of a certain equation, 
which describes the aircraft's short-period motion, fixed. This is done by 
using feedbacks with variable gains. The gain-changing mechanism is found 
using gradient techniques. 

The system was shown to be practical by an analog simulation. It was 
found to be tolerant of instrument noise, elevator hysteresis, and other com- 
plications not accounted in the analytical derivations. 

The study strongly suggests a modification of the C*-Criterion. It is 
proposed that the coefficients in the expression of the quantity C* should 
be representative of the aircraft's parameters rather than have the values 
currently in use. 
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FOREWORD 



Three Master's theses, currently In progress, are extensions of the 
studies outlined in this report and will use it as reference. The theses are: 

"An Investigation of Oscillations in an Adaptive Aircraft Control System under 
Large Input Commands," by Ens. L. S. Wisler, USN. 

"A Study of an Adaptive Aircraft Control System in a Self-Organizing Con- 
figuration," by Lt. H. M. Richarde, USN. 

"The Extension of an Adaptive Aircraft Control Concept to Helicopters, " by 
Lt. J. M. Hood, USN. 

Most of the work reported was done while the author was employed by 
Honeywell, Inc. The analytical derivations are presented but only the con- 
clusions derived from supporting studies with an analog computer are given. 
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Adaptive Control for Aircraft 



The usual method for accommodating an automatic flight control system to 
the wide variations in dynamic characteristics of the airframe, with changes of 
airspeed and altitude over its flight envelope, is to change control-loop gains 
with measurements of air-data. The dynamic pressure or the Mach number, as 
estimated by an air-data computer, is generally used as the parameter with 
which the gains are scheduled. The electronic art, developing at a bewilder- 
ing pace, is used to provide reliable, high-performance control hardware. All 
flight control systems currently in production, with the exception of an adap- 
tive system for the F-lll and some fixed-gain controllers, use air-data 
scheduling. This includes the system for the giant C-5A and flight controls 
proposed for the new cycle of commercial transport aircraft. However, the 
engineering analysis required for design is extensive and the final adjustments 
to tailor the control system to the aircraft during flight tests are expensive, so 
there has been the hope that methods of control would be found which did not 
rely on scheduling with air-data but changed parameters by measurements of 
dynamic performance. These would ultimately lead to much simpler, more reli- 
able and more universally applicable systems. 

A large number of such schemes has been proposed in the last fifteen 
years. Some have been studied in flight tests and only one has gone into pro- 
duction. A brief survey of the older ideas is given by Blakelock, Reference 1. 
Boskovich and Kaufmann, Reference 2, describe the evolution of a successful 
high-gain, wide bandwidth, model-follower system. The Bendix system, 
described by J. Bernadyn, et al. , Reference 3, was flown. The North American 
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SIDAC Controller, developed by Shipley and his associates. References 4 and 
5, has been extensively analyzed and completely designed but has been held 
short of flight test because of lack of funds. These devices, while being very 
successful in most respects, have debilitating characteristics. 

Controllers which detect limit cycles, either intentionally, as in the 
Honeywell device for the X-15, or unintentionally, as in the General Electric 
control system for the F-lll, are upset by air turbulence unless special pre- 
cautions are taken which then reduce the effectiveness and increase the com- 
plexity of the system. We conjecture that the Bendix OLAC was not carried 
further because it probably did not control adequately at all flight conditions. 
The SEDAC will probably show unwanted oscillations under large commands. 

It is a high-gain model-follower and is complicated. Since reliability, ease 
of maintenance, effect of component failures and results of battle damage are 
vital considerations besides performance in control system evaluation, the 
relatively complex adaptive schemes have been losing favor. Indeed, some 
people, for example Horowitz in Reference 6, argue that linear control should 
be adequate. 

A somewhat different point of view has been taken by R. C. Hendrick in 
developing a system, described in Reference 7, currently undergoing flight 
testing at Patuxent River. His approach is to use only the minimum amount of 
information from the dynamical motion that is needed for the single control sys- 
tem to be applicable to a whole class of aircraft without additional tailoring. 

Our study reverts to the goal of providing a response that is uniform for 
all flight conditions without using measurements of dynamic pressure, altitude 
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or angle-of-attack. Information obtained from inertial instruments, such as 
accelerometers or rate gyros, or from measurements of control surface posi- 
tions and command inputs is admissible. We persist in this area despite the 
general disenchantment with adaptive controls not only because it appears that 
a really practical solution for aircraft is close at hand but also because the 
study of the problem yields basic understanding of aircraft controls and the dif- 
ficult question of requirements for suitable handling qualities. 

The Approach 

Our scheme is very close to that used for SIDAC, (4, 5). The modifica- 
tion is based on the observation that the equation for handling-qualities C*- 
Criterion is very similar to a basic short-period equation for the motion of the 
aircraft. A modest feedback and feedforward with variable gains holds the co- 
efficients of the equation fixed. The C*-requirement may be met by choosing 
these coefficients to be the same as demanded by the Criterion or by adding 
a fixed outer loop. The mechanism for varying the gains is found by a gradi- 
ent calculation similar to the SIDAC-analysis . We achieve the advantage of 
a low-gain, narrow-bandwidth system which is very insensitive to instrument 
noise, bending modes, accommodates the primary control system and satisfies 
the requirements directly rather than using a model -following technique. 

The SIDAC-System identifies parameters and uses this in turn to adjust 
gains. Since the accuracy of identifying the several coefficients depends on 
the frequency content of the motion and the particular flight condition, it 
appears that a system which calls for the required response directly should 
have an advantage. This has also been argued by Hofmann and Best, Reference 
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8, in a fairly similar approach to control of the lateral-directional axes. 

The methods of synthesizing the gain changing mechanism have been 
developed by many authors. We studied a Lyapumov-function technique, re- 
ported in Reference 9. This was later rediscovered by Parks, Reference 10. 

P. M. Lion, Reference 11, has generalized and unified the several methods. 

The C*-Criterion 

The most difficult problem in flight control design, besides considerations 
of making the system invulnerable to component failures, is in deciding on what 
flying characteristics will be acceptable to the pilot. The requirement for longi- 
tudinal response which best fits into an analytical formulation is that promul- 
gated by Tobie, Elliot and Malcom, Reference 12. The criterion is that the 
time-response trace of a quantity called C*, for an abrupt force applied to the 
stick, must fall within a certain envelope, shown in Figure 1. C* is the sum 
of the normal force applied to the pilot's seat plus a constant multiple of the 
angular velocity in pitch. Thus, it is a linear combination, with constant posi- 
tive coefficients, of normal acceleration at the aircraft's center-of-gravity , the 
pitch velocity and the pitch acceleration, assuming the pilot's station is ahead 
of the eg. This blended response variable has been used for flight controls be- 
fore. It is a reasonable quantity to consider because at flight conditions with 
high dynamic pressures, at which the aircraft is very responsive, the pilot 
likes to feel the normal acceleration with stick force. He bases his assess- 
ment of the aircraft's performance on it. However, in low-q flight conditions, 
for example in landing, pitch rate is a more important clue to the aircraft's con- 
trollability. Indeed, for a helicopter at hover, pitch attitude is the variable 
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which must be controlled. The sum of the quantities exhibits this change in 



emphasis. The steady-state response of normal acceleration has the factor of 
the aircraft's velocity when compared to that for the pitch rate and hence nor- 
mal acceleration will dominate at high speeds and pitch rate will be important 
at low speeds. A more comprehensive analysis is given in the report by 
Makers, Reference 13. 

A Conjecture 

The study of our system by analog simulation showed that it was neces- 
sary to hold the coefficients of the short-period equation fixed as values 
representative of the ranges of the aircraft's coefficients over the flight enve- 
lope and not to achieve the C*-values directly. Analysis presented later also 
shows this. It is necessary because of higher order effects that are ignored 
in the basic derivation but are present in realistic situations. The C*-criterion 
is met by using an outer-loop feedback with fixed coefficients . 

This naturally prompts a review of the C*-requirement. For the F-4 air- 
craft, which was used as the vehicle for the study, the pilot sits only twelve 
feet ahead of the center-of-gravity . This makes the effect of pitch accelera- 
tion on the C*-quantity very weak. Since most of the studies of C* were done 
on fixed-base simulators, it appears possible that the angular-acceleration 
clue is not properly accounted in the criterion. Some studies, for example, 
done by Bihrle, Reference 14, indicate that it is an important clue. Thus, we 
conjecture that a more natural criterion for longitudinal handling-qualities is 
that the coefficients of the C*-expression be taken in the range of the coeffi- 
cients of the corresponding short-period equation. 
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Not only does this make our adaptive scheme more natural but it appears 



that the envelope of allowable responses. Figure 1, may be made much more 
restrictive. Then the Criterion becomes more definitive and easier to apply. 
Further, it can be directly extended to include helicopters by adding a term 
for pitch attitude. The corresponding short- period equation for the helicopter 
in hover has just such a term and this quantity is an important handling cue in 
this condition. Thus, this modification of C* is intuitively appealing. 



We begin with the equations of perturbations from straight and level flight, 
written as 



in terms of angular rate in pitch £, angle-of-attack of , and normal accelera- 
tion at the center-of-gravity ft. The quantity is the value of the aircraft’s 
unperturbed velocity and represents the elevator deflection. The coeffi- 



and have appropriate dimensions. All angles are measured in radians and 
linear quantities are measured in feet. Data for the F-4 aircraft is given in 



Equations for Short- Period Motion 




6 * M- « ♦ 

Y\ « Lfe (©-of) > 



(l) 




Table 1. 
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Angle-of-attack is difficult to measure, and its interpretation is compli- 
cated by turbulence in the air, so, following Shipley (4), we eliminate it from 
the equations and find 



e - p^©+pv,vi+p s s 

u o 



( 2 ) 



The coefficients are 

C js,= + M^. 
) - m V 2 « - 

ft • M s - *« 



and are listed in Table 1. They are negative for all flight conditions. The 
first equation in Set 2 is the fundamental relation in this study. 

The C*-quantity is 

C* - + 4 6 + ^ (4) 

in which £ is the distance of the pilot forward of the eg and U c is a number 
called the cross-over velocity, usually taken around four hundred feet per sec- 
ond. If we require C* to be precisely a multiple of the command input C c , 
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Equation 4 may be written as 




(6) 



This is exactly the form of the fundamental equation in (2). 

The Basic Derivation 

The control configuration is diagrammed in Figure 2. To avoid complicat- 
ing the discussion, for the present neglect the actuator, actuator model, the 
fixed feedback and the dynamics of the primary control system. Thus, assume 
for now that 




s 

c 



c ♦ -f 

C - C c 4 



(7) 



We are looking for a control function 



i 



in which the gains 






+ a, vi * v t c 

vary so that for each flight condition 



(8) 



the over -all system behaves as if the aircraft equation 



= o 



Q — ' V) - % rs 

were replaced by one with given fixed coefficients ft- Pi 



( 9 ) 






© “ ^ d “ ^V| ^ ~ O* 



( 10 ) 



Hence, we define an error signal, on which the generation of the variable 
gains is based, by the amount with which Equation 10 fails, namely 



0 — ^ © ** -> 



Subtracting (9) from (11) and incorporating the control function yields 



z « + §$ r^. j © * ^ p* p* 

+ c • 



v? 



(id 



(12) 



Therefore, to hold the error zero the gains must be adjusted to the values 



p* _ Syfr 

' fc 



i **% 



C*- , ff- (13) 

Pi 8 ' 
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Values of these with 

-3 , p^a-.oa. p, = -zo 

are listed in Table 1. 

To construct the mechanism for calculating the gain, we follow Shipley in 
considering the gradient of the sum of squares of the coefficients of Equation (12). 
This method is now a standard technique. Let 



where 




C + P* 

< ^ ^ + & n, 



(14) 



(15) 



and . VC,, ^ 



are positive numbers. 



Then differentiation of (14) with 



times gives 



av 

M 




4£n 

a-t 






<h j. 



(16) 



Hence, choosing 



arv. * 

if - 

a . - 

di ~ 

—nr - c ^ 



(17) 
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results in the equation 







(18) 



Since , choosing £ as 




c ® s 



(19) 



or 



r + » , 

Q * Sjv, t * V o , -fo i £ * £> 

k. - 1 , £ < - t» 



(2 0) 




the gains will tend to the ideal values at each flight condition. The system is 
not stable in any classical sense because the manner in which the gains con- 
verge depends on the activity and frequency content of the input. However, a 
concept of eventual stability has been used by LaSalle and Rath, Reference 15, 
to describe the mathematical behavior of such systems. 

We prefer the mechanization given by (2 0). Shipley at first used (19) and 
found that an automatic gain-control circuit was needed to provide stability 
under large inputs. The final SIDAC system (5) uses not only the sign of the 
error but also the sign of the dependent variables to form the quantities which 
when integrated produce the varying gains. In our system the gains are found 
from the equations 
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= VC* ^ v> %<* £ el-t 



k. ^ © S$Y1 £ «u 



( 21 ) 



with 




* k, C c Sj* £ ^ 

kv, * s«» f , k s ■= i 



Soo j €* = -°Z' 



Modifications to Include the 
Dynamics of the Actuator 

Returning to Figure 2 we now include the actuator, considered as a first 
order lag, add the corresponding model of the actuator, and the fixed-feedback 
outer-loop. The equations are then modified to read 



Since the Laplace operator is equivalent to differentiation, the first two equa- 
tions may be rewritten as 




(22) 



C ® C c ^ © 4* $ + Hh ^ * 



CT^s»oe * (tu>o® -p„r>-p»c 

(liS+t) ® s "It ^ d t 0 + T* pi'M > p»iVi 



(23) 



+ P» (c+-f )• 
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Combining these two and including the evaluation of the control function (8) 



gives 

(TaS+i) £ — T«.(p^- + 

*• *- fcfipe ♦ (Pn-#H ♦ fcoO" (24) 

+ ( Ps - P* **■ P$ ^5 ) c . 



Thus, the actuator introduces the first two terms on the right hand side which 

were not found in Equation 12. These could be taken into account by adding 

«# 

two more feedbacks of & and f> with varying gains. We avoid this compli- 
cation and hope that these terms will not lead to serious errors. Analog simu- 
lation shows that this is the case. 

The requirement given by Equation 5, that C* should be a multiple of the 
command input is unrealistic, considering the vast latitude by Figure 1. We 
modify this by requiring that C* follow the command by a simple lag. 




T c S + \ 




(25) 



and calculate the deviation from this by 



?c = c 



S -M 



Cc . 



(26) 



This is the equivalent to 

(TcStOf c = "Tc V) + V) + 4 Tc S + W+UcTe)^ ( 27 ) 

+ f |<C c 
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when the definition of C*, Equation 4, is used. The first equation of (23), upon 
adding the fixed feedback, becomes 



(iVS+OC - T+P„y> + 

" (Pi + P*Hv- ) v > - P» C c • 



(28) 



(29) 



Multiplying Equation 28 by a constant and subtracting it from Equation 27 
yields 

(TcS+l)£c as (Tc 

«+• (Tc ♦yw-Ti. v> 4-C^Vl^Uc -/V* 1 * P* 

+ ( <Jfc +/* (3f H* )0 *♦* ( I V* P* ^ 

*t* ( k.+^X< f5% ) Cc . 

We now have the quantities , W<* , y*. and k. at our disposal to make 

the C*-errar small. It is assumed that the gain changers hold the performance 
error £ small. Thus, require 

V 4 f 5 * ® ° 

I + +• - O 

Uc + yU 4- ^ PS ® 0 

JL + Tc Uc ~^y u - f / T ' “*>* Ps Mj = o . 

The first equation gives us a gain relation, the others define the fixed-feedback 



(30) 



gains as 



C H„ = - 



= - 



U pv. 



Uc 4- >< Pv 






“ ~ 



ps 

J?+7c Uc Py 

><Pi 



(31) 
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The quantity^ can now be taken as a gain for a root locus study or to kill off 
one of the two remaining terms in Equation (29). With the data 

1 = \Z To. = O.l = -20 > 

Tc — .IS" — 3 

u c - 3oo = -.Oi 

the variations of these gains, the open-loop steady-state gain G , and the re- 

o 

suiting closed-loop roots are shown in Table 2, (the ideal feedbacks (13) are 
assumed). We note that the open-loop gain is not very large. Bode plots show 
that the bandwidth is always less than 12 radians per second. This analysis 
was made from the transfer functions recorded in the Appendix. 

Performance of the System 

The controller was evaluated on an all-analog simulation. After initial 
study with only the actuator as the complication, a representation of the pri- 
mary control system, a second-order, forty-radian-per-second servo, a simula- 
tion of 0. 1-degree elevator hysteresis, and equations describing a bending 
mode were added. The system accepted these although the rates and accura- 
cies for the convergence of the gains were reduced. The dynamics of the 
second-order servo introduced limit cycles under relatively large command in- 
puts. The effects of noisy measurements and air turbulence were studied. The 
system was not upset by these. The most notable aspect of the performance 
was the insensitivity of the C*-response to large deviations of the varying 
gains from their ideal values. 

The fixed outer feedback could be eliminated by choosing the parameters 
of the error expression (11), f?* , ^ ^ to be equal to the corresponding 
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values in the C*-equation (16). The simulation showed that with this choice 



the gain-changing mechanism was unstable. To keep the loop-gains and the 
errors introduced by neglected complications small, these parameters are picked 
to be in the range of the corresponding aircraft parameters ' P*» » P* • 

Then the system is stable. 



Conclusion 



The Nonvarying-C* Control Scheme has been shown to have practical appli- 
cation for fixed-wing aircraft. It should be possible to extend the method for 
use on helicopters. The approach of choosing the variable-gain mechanism on 
the basis of system performance rather than requiring an explicit parameter 
identification results in a low-bandwidth system that is very tolerant to noise 
and high-order effects. The study also suggests a modification of the C*- 
Criterion for longitudinal handling qualities. 
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FIGURE 2. THE NONVARYING -C* CONFIGURATION 



Flight 

Cotop. 


1 


z 


3 


A 


S' 


4 


7 


6 


A|i. xio 3 


5L 


IS 


3o 


sl 


SC 


IS 


3o 


SO 


MocU Mo. 


0t2 


0.4 


as 


a 9 


1.2 


I.S 


2.0 


2.0 


Uo 


223 


-423 


591 


looS 


1340 


IS 87 


1990 


1938 


- M* 


1.3 


2.93 


3.4 


38 


180 


IS 7 


no 


42.4 


- 


*26 


.33 


.27 


1.38 


1.33 


.40 


. 10 


.04 






.54 


.si 


24| 


3.84 


240 


1.21 


<48 


” Ms 


z.a 


4.3 


7.5" 


SB. 4 


loi 


45. 0 


38 4 


14.9 


-2< 


83 


2Z8 


242 


2420 


4421 


3KZ 


2119 


817 




13. 2 


3I.S 


38.2 


313 


400 


401 


24 S“ 


102 




2. SB 


S'9o 


4.?J 


57. S 


74 4 


45/ 


24 8 


941 


-h 


,4? 




.7? 


AO 


5 2 


2-8 


13 


.52 


- ^ 


.01+4 


.OI2J 


.0133 


.0143 


• OJ97 


.<*94 


.0577 


.os/ 7 


r * 


.40 


.34 


.32 


- oi<r 


-.Ozd 


,oo4 


.048 


.24 




.0021 


. 00 13 


.00 10 


.0001 


-.0003 


-.0 004 


-.0013 


-.0033 


r s* 


6 . 75 * 


2.37 


/.88 


-.43 


-.74 


-.53* 




/.08 



TABLE 1. AIRCRAFT PARAMETERS 
(F-b Aircraft. Units in radians, feet, seconds.) 
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ROOT AND GAIN VARIATIONS WITH 



APPENDIX 



Transfer Functions 



Assume the ideal gains (13) obtain. Then the open- loop transfer function, 
breaking the loop before the actuator, has the numerator 

No - ^ C6fc Mv, + 

- + uo ( e s h* + pv, )] 

with denominator of the aircraft and the actuator as 

= (Ta.41-0 (fry 1-^ ) s + z* (($»,+ ^ )] ') 

the open- loop steady- state gain is 

r - + )+ U> (.PsH* 

Gto ~ ~ + <->„ 

Closing the loop produces the denominator 

A c = T a S S + 4 ^ )- pj(l- ) W^. + |r p-ojs’’ 

+ tSv"* ■ KT * a,,(p, ’ + ^ )+ 

- (* H « + P*>-P«d S +- |* [Mt + p t +u.((StH M + pj\. 



After a little algebra, the C*-transfer function is calculated to be 

C* _ ^ fr _ (E.SN- E.S + E. )S 

C c Tc S + 1 ^ c 



with 



Ea _ Uo CTc 1 ->*Ta (3v, )J 

E, =y<.T*. (s-T*-Tcjy+ 12 ft J s 

E t = H J? - Te ^ (' + -^Pv, ) - 
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